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Abstract. We have performed hydrodynamic simulations of relativistic rotational supernova core collapse in ax- 
isymmetry and have computed the gravitational radiation emitted by such an event. The Einstein equations are 
formulated using the conformally flat metric approximation, and the corresponding hydrodynamic equations are 
written as a first-order flux-conservative hyperbolic system. Details of the methodology and of the numerical code 
have been given in an accompanying paper. We have simulated the evolution of 26 models in both Newtonian 
and relativistic gravity. The initial configurations are differentially rotating relativistic 4/3-polytropes in equi- 
librium which have a central density of 10^" g cm~^. Collapse is initiated by decreasing the adiabatic index to 
some prescribed fixed value. The equation of state consists of a polytropic and a thermal part for a more real- 
istic treatment of shock waves. Any microphysics like electron capture and neutrino transport is neglected. Our 
simulations show that the three different types of rotational supernova core collapse and gravitational waveforms 
identified in previous Newtonian simulations (regular collapse, multiple bounce collapse, and rapid collapse) are 
also present in relativistic gravity. However, rotational core collapse with multiple bounces is only possible in a 
much narrower parameter range in relativistic gravity. The relativistic models cover almost the same range of 
gravitational wave amplitudes (4 x 10"^^ < h^'^ < 3 x 10~^° for a source at a distance of 10 kpc) and frequencies 
(60 Hz < f < 1000 Hz) as the corresponding Newtonian ones. Averaged over all models, the total energy radiated 
in form of gravitational waves is 8.2 x 10~* Mqc? in the relativistic case, and 3.6 x 10~* Mqc? in the Newtonian 
case. For all collapse models which are of the same type in both Newtonian and relativistic gravity, the gravita- 
tional wave signal is of lower amplitude. If the collapse type changes, either weaker or stronger signals are found 
in the relativistic case. For a given model, relativistic gravity can cause a large increase of the characteristic signal 
frequency of up to a factor of five, which may have important consequences for the signal detection. Our study 
implies that the prospects for detection of gravitational wave signals from axisymmetric supernova rotational core 
collapse do not improve when taking into account relativistic gravity. The gravitational wave signals obtained 
in our study are within the sensitivity range of the first generation laser interferometer detectors if the source 
is located within the Local Group. An online cat alogue containing the gravitational wave signal amplitu des and 



spectra of all our models is available at the URL http : //www. mpa-garching .mpg . de/Hydro/hydro .html 
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1. Introduction 

At the end of their thermonuclear evolution, massive stars 
develop a core composed of iron group nuclei (hence iron 
core) which becomes dynamically unstable against grav- 
itational collapse. The iron core collapses to a neutron 
star or a black hole releasing gravitational binding energy 
of the order - 3 x IQ^'^ erg {M/Mq)'^{R/ 10 km)-\ which 
is sufficient to power a supernova explosion. If the core 
collapse and/or the supernova explosion are nonspherical, 



part of the liberated gravitational binding energy will be 
emitted in the form of gravitational waves. Nonsphericity 
can be caused by the effects of rotation, convection and 
anisotropic neutrino emission leading either to a large- 
scale deviation from spherical symmetry or to small-scale 
statistical mass-energy fluctuations (for a review, see e.g., 
MiiUer (1998|)). 



According to present knowledge the energy radiated 
away in form of gravitational waves does not exceed 
IO^^A/qC^ in rotational core collapse (MiiUer 1982; Finn & 
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Evans 1990; Monchmcyer et al. 1991 



Bonazzola & Marck 



1993; Yamada fc Sato 1994a; Zwerger fc MiiUer 1997 



2 



H. Dimmelmeier, J. A. Font & E. Miiller: Relativistic rotational core collapse. II 



Rampp et al. 1998; Fryer et al. 2002) and from convection 
or anisotropic neutrino emission in neutrino-driven super- 
novae ( Miiller fc Janka 1997] ) . The frequency of the emit- 
ted radiation ranges from about a few Hz to a few kHz, 
and the (dimensionless) signal amphtudes for a source lo- 
cated at a distance of 10 Mpc (within the Virgo cluster) do 
not exceed ~ 10~^^. The smallness of the signals calls for 
accurate waveform templates to alleviate the technical dif- 
ficultie s in extracting the signal from noisy data ([Pradier 



et al. zOOl). For the current generation of laser interfero- 



Here il, I, Epot, Rs, and R are the (constant) angular 
velocity, the moment of inertia about the center of mass, 
the gravitational potential energy, the Schwarzschild ra- 
dius, and the radius of the star, respectively. The positive 
constant k depends on the density distribution of the star. 

Axisymmetric Newtonian hydrodynamic collapse and 
explosion simulations using a realistic equation of state 
(EoS) and some treatment of weak interaction processes 
have been performed by Miiller (1982| ) 



Bodcnheimei 



metric detectors (GEO 600, LIGO, VIRGO, TAMA) such 



small amplitudes imply that the prospects for detection 



of gra\ itational waves from core collapse supernovae are 



fc Woosley (1983) ), |Symbahsty (1984| ) and [Bonazzola 
fc Marck (1993|) neglecting ne utrino transport, and by 
Monchmeyer fc Miiller (1989|) , Janka & Monchmcyei 



limited to those events occurring within the Local Group. 
However, if measured, a gravitational wave signal does 
provid e (as observations of neutrinos) a direct diagnosis 



Nadezhin 



(1989[) ,' |Monchmeyer et al. (1991D , llmshennik 

(T992I , pV^r fc Hegcr (20001 ) and pryer et al. (2002| ) cm- 



ploying some approximative description of neutrino trans- 
port. In addition. 



of the lynamics of the events. 



Finn fc Evans (1990|) , [Yamada fc Sate 



(19941 ) and ^werger fc Miiller (1997 ) have performed 



Investigations of rotational core collapse are also im- 
port a nt in the context of the supernova explosion mecha- 



Janka fc Monchmeyer 1989; Imshennik fc Nadezhin 1992; 



1983|; ^ymbahsty 1984|; [Monchmeyer fc Miiller 1989 



Newtonian parameter studies of the axisymmetric collapse 
of rot ating polytropes. Rampp et al. (199^) and [Brown 



nism (| /liillcr fc Hillcbrandt 1981 Bodenhcimer fc Woosley (200l|) have extended the work of [Zwerger fc Miiller (1997| ) 



by relaxing the assumption of axisymmetry. 



Yamada fc Sato 19941] ; [Fryer fc Hcger 2000[ ) and for the 
coUapsar scenario of gamma-ray bursts, where the collapse 
of a rotating massive star leads to the formation of a Kerr 



black hole and a relativistic jet (MacFadyen et al. 2001 
Aloy et al. 2000] ; ]Wheeler et al. 2000^ 



Wilson (1979), [Evans (1986[ ), [Nakamura (198l[ , [1983] ) 
and Stark fc Piran (1985[ ) pioneered investigations of 
axisymmetric collapse of rotating configurations in full 
general relativity (GR). Wilson (1979[ ) computed neu- 
tron star bounces of 7 = 2 polytropes, while Nakamura 
(1981] ) (see also Nakamura et al. (1987)) simulated the 



Numerous studies have addressed quite different as- 
pects of rotational core collapse at different levels of so- 
phistication. Ideally, such studies should incorporate a rel- 
ativistic treatment of gravity and hydrodynamics with- 
out symmetry restrictions, a detailed description of the 
complex microphysics and of the neutrino transport, and 
should use consistent initial models from evolutionary cal- 
culations of rotating stars. The latter point is still a nag- 
ging problem (for a major step forward in this context see 
Heger et al. (2000[)), as up to now all studies, except that 
of Fryer et al. (200^ ) , have relied on parameterized initial 
models. 

As the general relativistic corrections of the gravita- 
tional potential during supernova core collapse, and for 
a neutron star, do not exceed 30%, it is often argued 
that they can be neglected or only taken into account 
approximately. While this approach may be justified to 
some degree in the case of non-rotational core collapse, it 
becomes very questionable when the core possesses a sig- 
nificant amount of angular momentum. In the latter case, 
a relativistic treatment of gravity is much more impor- 
tant, because the stabilizing effect of rotation is counter- 
acted by the destabilizing effect of the deeper relativistic 
potential. This is refiected in the expression for the crit- 
ical adiabatic index 7crit, below which a rigidly rotating 
relativistic configuration is dynamically unstable against 
pseudo-radial (linear) isentropic perturbations (see, e.g.. 
Chap. 14 of [Tassoul (l"978[ )): 



formation of rotating black holes resulting from the col 
lapse of a IOMq "core" of a massive star with different 
amounts of rotational energy and an initial central density 
of 3 X 10^'^ g cm"'^. Nakamura (1983) (see also Nakamura 
]et al. (1987 )) considered a configuration consisting of a 
neutron star (M — 1.O9M0, pc — 10^^ g cm^^) with an 
accreted envelope of O.8IM0, which was thought to mimic 
mass fall-back in a supernova explosion. To this configu- 
ration he added rotation and infall velocity, and simu- 
lated the evolution depending on the prescribed rotation 
rates and rotation laws. In both scenarios the EoS con- 
sisted of a relativistic degenerate lepton gas (7 = 4/3) at 
low densities {p < p* = 3 x 10^^ g cm~^, and of a stiff 



(7 = 2) com ponent at large densities {p > p*). Stark fc 



Piran (1985 ) were the first to compute the gravitational 
radiation from the relativistic collapse of a rotating poly- 
tropic (7 — 2) star to a black hole. The initial model was a 
spherically symmetric relativistic polytrope in equilibrium 
of mass M, central density 1.9 x 10^^(M/Mo)~^, and ra- 
dius GGM/c^ = 8.8 X 10^M/Mq cm. Rotational collapse 
was induced by lowering the pressure in the initial model 
by a prescribed fraction, and by simultaneously adding 
an angular momentum distribution approximating rigid- 
body rotation. Stark fc Piran (1985) found a low efficiency 
of gravitational wave emission (-Erad tot M^c^ < 7 x 10~'*, 
where E'rad tot is the total energy radiated by graviational 
waves), and that for sufficient rotation the star bounces 
and no black hole forms. 



_ 4 2 n'^I Rs 
- 3"9^ + ^l^■ 



(l) 



After the work of Stark fc Piran (1985) it took about 
15 years before the next simulations of GR rotational core 
collapse were published. This was mainly caused by per- 
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sistent numerical problems occurring in axisymmetric GR 
simulations due to coordinate singularities at the symme- 
try (= rotation) axis. These coordinate singularities ham- 
per the development of methods which guarantee accuracy 
and stability in long term (covering many dynamical time 
scales) simulations. Alcubierre et al. (2000| ) proposed a 
method which does not suffer from stability problems and 
where, in essence, Cartesian coordinates are used even for 



axisymmetric systems. Using this method Shibata (2000) 
investigated the effects of rotation on the criterion for 
prompt adiabatic collapse of rigidly and differentially ro- 
tating (7 = 2) polytropes to a black hole. Collapse of the 
initial approximate (computed by assuming a conformally 
flat spatial metric) equilibrium models was induced by a 
pressure reduction. [Shibata (2000 ) found that the criterion 
for black hole formation depends strongly on the amount 
of angular momentum, but only weakly on its (initial) 
distribution. He also studied the effects of shock heating 
using a gamma-law EoS, and found that shock heating is 
important in preventing prompt collapse to black holes in 
case of large rotation rates. 

Hayashi et al. (1999) investigated the possibility of sec- 
ular, i.e. quasi-static, core contraction from white dwarf 
to neutron star densities using equilibrium sequences of 
rapidly rotating, general relativistic compact stars with 
phcnomcnological equations of state. They demonstrated 
that there is a possibility for the existence of "fizzlers" in 
the framework of GR, at least for the simplified equations 
of state used in their study. 

The above discussion shows that previous investiga- 
tions in GR rotational core collapse were mainly con- 
cerned with the question of black hole formation under 
idealized conditions, but none of these studies has really 
addressed the problem of supernova core collapse which 
proceeds from white dwarf densities to neutron star den- 
sities, involves core bounce, shock formation, and shock 
propagation. Exactly this is the motivation for the present 
work. To this end, we apply the numerical methodology 



presented in an accompanying paper (Dimmelmeier et al 



2002, hereafter Paper I) and study the dynamics of ax- 
isymmetric, relativistic rotational supernova core collapse 
and the associated emission of gravitational radiation. In 
particular, we have simulated the collapse of a compre- 
hensive set of rotating relativistic stellar models param- 
eterized by the initial degree of differential rotation, the 
initial rotation rate and the adiabatic indices at subnu- 
clear (supranuclear) matter densities. Thus, our study is 
a GR extension of the previous work of ^werger (1995 ) 
and Zwerger & Miiller (1997), who computed the gravita- 
tional radiation from rotational core collapse supernovae 
using Newtonian gravity. First results from our study 
have already been published in a short communication 
(Dimmelmeier et al. 2001)]. 



The organization of the paper is as follows: Section || 
contains a brief repetition of the assumptions made in our 
study (for more details see Paper I). In Section ^ we an- 
alyze in detail the collapse dynamics and introduce the 
three different types of collapse and gravitational wave- 
forms we have identified. Section^ is devoted to the grav- 
itational wave emission of the models, and to the prospects 
of detectability of gravitational waves from core collapse 
supernovae. We conclude with a summary in Section ^ 
Additional information is given in Appendix ^ where we 
discuss the issue of gravitational wave extraction. 



2. Model assumptions 

The matter model obeys an ideal gas EoS with the pres- 
sure P consisting of a polytropic and a thermal part for 
a more realistic tre a tment of the effects of s hock waves 
( [lanka et al. 1993| ; ^werger fc Miiller 199"?! ). Since we 



^ Note that due to a programming error, which we found only 
after the Letter was published, the values of some quantities 
(e.g., the bounce densities and maximum wave amplitudes) 



are mostly interested in the gravitational radiation emis- 
sion, which is controlled by the bulk motion of the fluid, 
we neglect microphysics like electron capture and neu- 
trino transport. The initial conflgurations are differen- 
tially rotating relativistic 4/3-polytropes in equilibrium 
( Komatsu et al. 1989a| ,[b| ^tergioulas fc Friedman 1995 ) 
which are marginally stable, and which have a central den- 
sity pcini — 10^*^ g cm~^. Collapse is initiated by decreas- 
ing the adiabatic index (initially 4/3 in all models) to some 
prescribed flxed value 71 (for more details see Paper I). 

The initial models are determined by three parame- 
ters, which are also used to name the models. The first 
parameter A is a length scale, which specifies the degree of 
differential rotation (see Eq. (33) of Paper I). The smaller 
the value of A the more differentially rotating the model 
is. Model parameters Al, A2, A3 and A4 correspond to 
A = 5 X 10^ cm, A = 10^ cm, A = 5 X lO"^ cm, and 
A = 10^ cm, respectively. The second parameter is the 
initial rotation rate /3rotini, which is given by the ratio of 
rotational energy and the absolute value of the gravita- 
tional binding energy. Model parameters Bl, B2, B3, B4 
and B5 correspond to /3rotini = 0.25%, 0.5%, 0.9%, 1.8% 
and 4%, respectively. The third parameter 71 is the adia- 
batic index at subnuclear densities (p < Pnuc), with model 
parameters Gl, G2, G3, G4 and G5 corresponding to 
71 = 1.325, 1.320, 1.310, 1.300 and 1.280, respectively. The 
name of a simulated model is then given by a combination 
of parameters from the three sets (e.g., A3B4G5). We de- 
fine the nuclear matter density p^uc as 2.0 x 10^** g cm~^. 

The adiabatic index at supranuclear densities (p > 
Pnuc) is fixed to 72 = 2.5, except in one model which has 
also been run with 72 = 2.0 in order to test the influence 
of a softer supranuclear EoS on the collapse dynamics. 

In total, we have simulated the evolution of 26 models 
(see Table 0). In order to identify the relativistic effects 
on the collapse dynamics, we have simulated these models 
also in Newtonian gravity. Note that the Newtonian simu- 



given in the Letter differ from those presented here, but all 
qualitative results remain unchanged. 
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lations are a subset of those performed by Zwerger (1995) 
and [Zwerger fc Miiller (1997| ) . 

The Einstein equations are formulated using the 
so-called conformally fiat (CF) metric approximation 
( Wilson et al. 1996 , conformal flatness condition - CFC) 
and the corresponding hydrodynamic equations are formu- 
lated as a first-order flux-conservative hyperboUc system 
(Banyuls et al. 1997), well-adapted to numerical schemes 
based on Riemann solvers. The applicability and quality 
of the CFC for rotational core collapse has been discussed 
in detail in Paper I, where we demonstrated that its usage 
is appropriate. 



3. Collapse dynamics 

3.1. Collapse and waveform types 

Zwerger (1995) first identified three distinct classes of core 
collapse types. This classification is based on the form 
of the gravitational wave signal. However, as the gravita- 
tional radiation waveform is closely linked to the collapse 
dynamics, it also mirrors the collapse behavior. Thus, the 
signal types can generally be used to classify the collapse 
type, too. 

In our relativistic simulations, each of the examined 
models belongs to one of Zwerger's three collapse types, 
which he called type I, II and 111. In the following we 
describe a representative model from each class, discuss its 
characteristic properties and explain the physical effects 
which lead to its collapse classification. 



settles down into a new equilibrium state, which we will 
call the "compact remnant"^ in the following. While the 
shock wave propagates outwards, the inner core oscillates 
for about 10 ms with a superposition of several damped 
eigenmodes with frequencies of about 10'^ Hz (ring-down). 
If the infall is moderately fast (collapse timescale longer 
than about 50 ms), we will talk of a regular collapse model 
giving rise to a type 1 gravitational wave signal according 
to the nomenclature of Zwerger & Miiller (1997). In case 
of a very rapid infall (collapse timescale of about 30 ms), 
the core plunges unimpeded by centrifugal forces deeply 
into the gravitational potential well, with its central den- 
sity significantly exceeding nuclear density. In such rapid 
collapse models, the gravitatio nal waveform (type 111 ac - 
cording to the nomenclature of Zwerger fc Miiller (1997 )) 
is qualitatively different from that of a regular (type I) 
coUapse model (see below). 

On the other hand, if core collapse is only or predom- 
inantly stopped by centrifugal forces, the inner core ex- 
periences several distinct sequences of infall, bounce, and 
re-expansion separated by up to 50 ms before it eventu- 
ally settles down into an equilibrium state. During each 
bounce an outward shock is generated. The multiple large 
scale bounces occur because the deceleration by centrifu- 
gal forces is less abrupt than that due to the stiffness of the 
supranuclear EoS. The latter is characterized by an adi- 
abatic index 72 = 2.5 (see Section ||), while rotation acts 
like a 7 = 5/3 gas according to the virial theorem (see, 
e.g., Tassoul (1978 ), Chapter 14). Multiple bo unce mod- 
els produce type 11 gravitational wave signals ( [Zwerger^ 



can be divided mto three phases (see, e.g., Muller (199fc JJT 
Infall phase: This phase covers the initial collapse of 
the core from the onset of the gravitational instability trig- 
gered by the sudden softening of the EoS due to the re- 
duction of the subnuclear adiabatic index. The inner part 
of the core, which collapses homologously (vr ''^ r), con- 
stitues the "inner core", while the "outer core" is falling 
supersonically. Depending on the model parameters the 
infall phase lasts between 30 ms and 100 ms. 

Bounce phase: When the EoS stiffens because repul- 
sive nuclear forces become important at densities above 
nuclear density pnuc = 2 x 10^** g cm~^, or when cen- 
trifugal forces begin to dominante over gravitational at- 
traction due to angular momentum conservation and ro- 
tational spin-up, the inner core decelerates on a timescale 
of about 1 ms. Because of its large inertia and infall kinetic 
energy, the inner core does not come to rest immediately. 
The core overshoots the equilibrium configuration it will 
eventually approach, and bounces back, which results in 
the formation of a shock wave at the outer edge of the 
inner core. During core bounce the amplitude of the grav- 
itational wave signal is largest. 

Ring-down phase, or re-expansion phase: If centrifu- 
gal forces remain sufficiently small until nuclear density 
is reached in the center of the core, the bounce occurs at 
central densities slightly in excess of nuclear density due to 
the stiffening of the EoS. In this case the inner core rapidly 



In general, the evolution of any supernova core collapse Miiller 1997D which consist of a distinct large amplitude 



peak for every bounce. 

If the peak central density at bounce is close to (or 
exceeds) nuclear matter density, and if the core rotates 
differentially and rapidly, a mixture of a regular and mul- 
tiple collapse type model can occur. In these transition 
models of type I /II the core re-expands less than in a gen- 
uine multiple bounce model, but it still exhibits distinct 
bounces and coherent large scale oscillations. 

3.1.1. Regular collapse - type I 



As in the work of Zwerger & Miiller (1997) we choose 
model A1B3G3 as a representative model for a regular col- 
lapse. Fig. 1^ shows the time evolution of the central den- 
sity pc and the gravitational wave signal amplitude Afg ; 
the definition of the wave amplitude Afp and the descrip- 
tion of the numerical wave extraction technique are given 
in Appendix ^ The three phases of the collapse are clearly 
identifiable in both panels. During the infall phase pc in- 
creases and exceeds nuclear matter density at 48.26 ms. 
The peak value of the density is reached at the time of 



^ Due to neutrino cooling and other microphysical processes, 
which we do not consider in our models, the compact remnant 
later evolves into a neutron star. This may be accompanied by 
mass accretion due to fall-back. We have not investigated this 
late time evolution of the compact remnant. 
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Fig. 1. Time evolution of the central density pc (upper 
panel) and the gravitational wave signal amplitude ylfp 
(lower panel) for the regular collapse model A1B3G3. The 
horizontal dotted line in the upper panel marks nuclear 
matter density /Onuc, and the vertical dotted line indicates 
the time of peak central density 
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Fig. 2. Time evolution of the central density pc (up- 
per panel) and the gravitational wave signal amplitude 
Afo (lower panel) for the multiple bounce collapse model 
A2B4G1. The vertical dotted lines mark the times of peak 
central density for each bounce. 



bounce at th — 48.63 ms with Pch = 4.23 x 10^^ g cm^^, 
which is about twice nuclear matter density. Subsequently, 
the core slightly re-expands and rings down to an equilib- 
rium state with a central density of pcf ~ l-Spnuc- 

The gravitational wave amplitude (lower panel of 
Fig. ^ increases during the infall phase, but has a negative 
peak value with |Afo |max — 659 cm at ig„ = 48.99 ms > 
t\j. In all regular collapse models the maximum signal 
amplitude is negative. For some models (not in model 
A1B3G3) the signal amplitude exhibits a small local min- 
imum around the time of bounce, a feature which was 



below nuclear matter density. Such subnuclear bounces 
can indeed be observed in some of our collapse models, 
e.g. model A2B4G1, for which the time evolution of the 
central density pc and the gravitational wave signal am- 
plitude A^Q are shown in Fig. ^. 

In both the central density evolution and the signal 
waveform distinct extrema are discernable. The (negative) 
maximum of the signal amplitude jAfp jmax — 548 cm is 
reached at a time igw = 101.16 ms, which is shortly after 



the time ty, = 99.78 ms when the central density reaches 
its first peak Pch — 0.90 x 10^"* g cm~'^ w 0.45pnuc- During 
the subsequent bounces the peaks of the central den- 
sity are significantly smaller than the peak value at first 
bounce. The prominent peaks in the gravitational wave 
signal are clearly associated with the distinct bounces vis- 
ible in the density evolution. 

find that during the ring-down phase the maxima of A^q We point out that model A1B3G1, which shows mul- 

are les s damped than the minima, because larger acceler- tiple bounces in Newtonian gravity (Zwerger & Miillei 



identified and discussed by Zwerger & Miiller (1997). The 
ring-down phase of the pulsating inner core is directly re- 
flected in the wave signal which oscillates accordingly. As 
in the Newtonian runs of Zwerger & Miiller (1997), we 



ations p,re encountered at the high density extrema of the 
core. 

3.1.2. Multiple bounce collapse - type II 

For an adiabatic index 71 close to the initial value 7ini = 
4/3 (see Section^, and for rapid and highly differential ro- 
tation, centrifugal forces can halt the collapse at densities 



1997 ), does not exhibit such a behavior in the relativistic 



simulations. This change of collapse type due to relativis- 
tic effects is discussed in Section 3.3. 



3.1.3. Rapid collapse - type III 

If the core collapses very rapidly due to values of 71 which 
are much smaller than 4/3, the qualitative behavior of 
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Fig. 3. Time evolution of the central density pc (upper 
panel) and the gravitational wave signal amplitude ylfp 
(lower panel) for the rapid collapse model A1B3G5. The 
horizontal dotted line in the upper panel marks nuclear 
matter density pnuc, and the vertical dotted line indicates 

the time nf peak rpntra.l density Note the qnalitativp 



diffpreirp in the wavpfnrm rnmparnd to the regular rnl- 

lapse (type I) model shown in Fig. |l|. 



the density evolution is very similar to that of the regular 
(type I) collapse model, except that the bounce and equi- 
librium densities are slightly higher, and the post-bounce 
pulsations of the inner core are even more strongly sup- 
pressed. For example, in model A1B3G5 (Fig. ||), which is 
a typical rapid collapse model, the values of the peak and 
the final central density are Pch = 4.55 x lO^"* g cm~'^ « 
2.5/3nuc (reached at = 30.20 ms), and pd ~ l.Spnuc- 

On the other hand, the shape of the gravitational wave 
signal amplitude is qualitatively different compared to 
type I models. This is shown in the lower panel of Fig. ^, 
where the time evolution of A^q is plotted. The maximum 



signal amplitude 1^20 | max = 124 cm at igw = 30.03 ms is 
now assigned to the first positive peak of the waveform 
contrary to the signal waveform of the regular collapse 
models. In the case of rapid collapse, the negative peak 
is almost entirely suppressed and the overall amplitude of 
the signal is significantly smaller. There is no clear fea- 
ture in the waveform allowing one to identify the time of 
bounce tb- 

For a given initial model, the transition between type I 
and type III collapse models occurs gradually with de- 
creasing values of 71 . In the gravitational wave signal this 




40.0 45.0 
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55.0 



Fig. 4. Time evolution of the central density pc (upper 
panel) and the gravitational wave signal amplitude Afo 
(lower panel) in the relativistic (solid lines) and Newtonian 
(dashed lines) simulation of the regular collapse model 
A3B2G4soft- The horizontal dotted line in the upper panel 
marks nuclear matter density pnuc • The vertical dotted line 
indicates the time of peak central density th- 



is reflected by a decrease of the negative main peak and 



an increase of the first positive peak (Zwerger & Miillei 



1997) 



3.1.4. Additional remarks 

In the above models the central density pc coincides with 
Pmax, the maximum density of the model at any time. 
However, very rapidly and differentially rotating mod- 
els, like all 'A4' models, either already possess a toroidal 
density distribution initially, or develop one during their 
evolution. In such models the maximum density is at- 
tained at some off-center location, and pc < Pmax holds. 
Nevertheless, the time evolution of Pmax can be used to 
classify the collapse type according to the criteria dis- 
cussed above. 

Our results show that relativistic simulations of ro- 
tational core collapse exhibit the same qualitative col- 
lapse dynamics and allow for the same classification as 



Newtonian simulations (Zwerger 1995; Zwerger & Miillei 



1997). However, there exist important differences in other 
aspects of the collapse arising from relativistic effects. 
These are discussed in the subsequent sections. 

3.2. Compactness of the core remnant 

For all models the central density at bounce pch (for all 
A4 models and model A3B5G4 the off-center maximum 
density Pmaxb) is larger in the relativistic simulation than 
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Table 1. Summary of important quantities for all models in relativistic (R) and Newtonian (N) gravity, t^, is the time 
of bounce, Pmaxb is the maximum density at bounce, lAfo |max is the maximum gravitational wave amplitude, fmax 
is the frequency of maximum spectral energy density, /?iotmax is the maximum rotation rate, Pmaxf is the maximum 
density after ring-down, and Type specifies the collapse type (for a definition see text). If the maximum density is 
located off-center, the respective density value is marked with an asterisk. For models of type II or I/II, Pmaxf cannot 
be determined. 
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Fig. 5. Same as Fig. ^, but for the multiple bounce col- 
lapse model A2B4G1. The vertical dotted lines mark the 
times of peak central density tb for each bounce. 



in the corresponding Newtonian one, the relative increase 
reaching ^ 700% in special cases (see Table 0). This also 
holds for the maximum density of the compact remnant 
after ring-down in case of instant formation of a stable 
equilibrium state (types I and III). This generic property 
is illustrated in the upper panel of Fig. ^ and Fig. || for 
the type I model A3B2G4soft (where relativistic effects are 
particularly large due to its soft EoS and thus high central 
density) and the type II model A2B4G1, respectively. 

When comparing models of the same collapse type in 
Newtonian and relativistic gravity, we find that although 
all relativistic models reach higher central and/or maxi- 
mum densities and, on average, larger infall and rotation 
velocities, only six models (see Table 0) have also larger 
maximum signal amplitudes than their Newtonian coun- 
terparts. The maximum signal amplitudes of all other rel- 
ativistic models are (up to 57%) smaller than those of 
the corresp onding Newtonian ones (T able |l|). As already 
discussed in Dimmelmeier et al. (2001 ) , the reduced maxi- 
mum signal strength can be explained by the fact that the 
amplitude, which is calculated using the quadrupole for- 
mula (see Appendix ^), is determined by the hulk motion 
of the core rather than just by the motion of the densest 
mass shells, strictly speaking by the second time deriva- 
tive of the quadrupole moment. Therefore, a core which 
is more condensed in the center can give rise to a smaller 
gravitational wave signal than a core which is less cen- 
trally condensed, but which is denser and moves faster in 
its outer regions. This is reflected by the weight factor 



in the integrand of the quadrupole formula ( A. 4 ) . Due to 
this factor matter at moderately large r contributes signif- 
icantly to the gravitational wave signal amplitude, while 
at very large radii (outside the inner core) the decrease in 
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Fig. 6. Radial density profiles Pc along the equator near 
bounce time for model A3B2G4soft (upper panel) and 
model A2B4G1 (lower panel). In both models, the den- 
sity of the relativistic simulation (solid line) is larger than 
that of the Newtonian simulation (dashed line) in the cen- 
tral region. At r « 6 km (upper panel) and r « 35 km 
(lower panel) the density of the relativistic model drops 
below that of the Newtonian one. The radii of the density 
crossing are indicated by the vertical dotted lines, and the 
horizontal dotted lines mark nuclear matter density Pnuc- 



density more than compensates the increase of the weight 
factor. 

The upper panel of Fig. ^ shows the equatorial density 
profiles in the inner core of model A3B2G4soft for both the 
relativistic (solid line) and the Newtonian (dashed line) 
simulation at bounce. As already seen in the central den- 
sity evolution (Fig. ^), the simulation in relativistic grav- 
ity yields up to 36% larger densities in the central regions. 
However, at radii r > 5 km the density in the Newtonian 
model exceeds that of the relativistic one. This "density 
crossing" , which takes place at different radii for different 
polar angles 9, can be observed in all the models we have 
investigated, at all times t > % and at all polar angles. 
The surface of density crossing is located inside the com- 
pact remnant for type I and type III models, and in the 
case of type II multiple bounce models in the coherently 
re-expanding and contracting inner core, which has not 
yet settled into a compact remnant at the end of the sim- 
ulation. The effect is persistent in time as demonstrated 
in Fig. 1^, which shows the time evolution of the density 
for model A3B2G4soft at two different radial locations in 
the equatorial plane. After bounce, the central density pc 
of the relativistic configuration is always larger than the 
corresponding Newtonian one (upper panel). However, at 
larger radii, but still inside the compact remnant (lower 
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Fig. 7. Time evolution of the density at different radii 
in the relativistic (solid lines) and Newtonian (dashed 
lines) simulation of model A3B2G4soft- Whereas the cen- 
tral density (upper panel) is larger in relativistic gravity, 
the (equatorial) density of the Newtonian model exceeds 
that of the relativistic model at a radius of r = 12.9 km 
(lower panel). The vertical dotted lines indicate the time 
of bounce tb, and the horizontal dotted line in the up- 
per plot indicates nuclear matter density pnuc- Note the 
different vertical scale of the two plots. 
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Fig. 8. Weighted radial density profile pcr^ along the 
equator (solid lines denote relativistic models; dashed lines 
denote Newtonian models) close to the time of core bounce 
for model A3B2G4soft (upper panel; t — 39.25 ms) and 
model A2B4G1 (lower panel; t — 99.75 ms), respectively. 
The crossing radii are indicated by the vertical dotted 
lines. 



panel), the density of the relativistic model is smaller by 
up to 50%. 

The analysis of the multiple bounce model A2B4G1, 
which yields a smaller signal amplitude in relativistic grav- 
ity (see lower panel of Fig. H) , leads to a similar result as 
the regular bounce model A3B2G4soft- In the relativis- 
tic case, the central density at bounce is more than 8 
times larger than in the Newtonian one. As in model 
A3B2G4soft, the equatorial density profiles cross (lower 
panel of Fig. |6|), but now the crossing occurs at a much 
larger radius r « 35 km. Hence, the region where the rel- 
ativistic model is denser is considerably more extended. 
Inside the crossing radius, the cores of both the relativistic 
and the Newtonian model oscillate coherently. The den- 
sity of the relativistic model is considerably higher in this 
region, and the matter experiences stronger accelerations, 
which causes narrower density peaks and shorter time in- 
tervals between two consecutive bounces (Fig. ^). 

The effects of the reduced density in the outer regions 
of the inner core in relativistic gravity on the gravitational 
wave signal can be inferred from Fig. ^, which shows the 
equatorial density profile weighted by a factor . Inside 
the crossing radius r « 6 km, the profiles of the weighted 
density pr^ of model A3B2G4soft (upper panel) coincide 
closely in both relativistic and Newtonian gravity. On the 
other hand, for 6 km ^ ^ 13 km the weighted den- 



sity pr^ of the Newtonian model is significantly larger 
than that of the relativistic model. Consequently, as the 
density crossing occurs at all polar angles, the resulting 



quadrupole moment (Eq. A. 4) is larger, too. In order to 
see whether the enhanced quadrupole moment also gives 
rise to a stronger gravitational wave signal, one has to con- 
sider the seco nd time derivative of the quadrupole moment 
(see Eq. A^), which results from the acceleration of the 
matter distribution in the core. We find that despite the 
somewhat lower average acceleration in Newtonian grav- 
ity which is reflected in the 18% lower frequency of the 
gravitational wave signal, the large quadrupole moment of 
the more extended density distribution results in a maxi- 
mum signal amplitude jvlfg |max, which is 26% larger than 
that of the centrally denser compact remnant in relativis- 
tic gravity (see lower panel of Fig. |^ and Table [l]). 

For the multiple bounce model A2B4G1, the situation 
is similar (lower panel) . In relativistic gravity the weighted 
density pr^ is enhanced in the central parts of the in- 
ner core, and reduced in the outer parts compared to the 
Newtonian profile. Contrary to model A3B2G4soft, the rel- 
ativistic density increase close to the center has a strong 
effect on the weighted density profile inside the crossing 
radius r « 35 km. Nevertheless, and despite the higher 
accelerations (i.e., larger second time derivatives) of the 
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coherently oscillating inner core in relativistic gravity, the 
smaller densities at larger radii again yield a smaller grav- 
itational wave signal (see lower panel of Fig. ||). However, 
in relativistic gravity the increase of both the central den- 
sity (718% vs. 36%) and the frequency (41% vs. 23%) are 
larger in model A3B2G4soft compared to model A2B4G1. 



1000.0 



Thus, the maximum signal amplitude 1^20 | max of the for- 
mer model is only reduced by 16% instead of 21% in model 
A2B4G1 when changing from Newtonian to relativistic 
gravity (see Table ||) . 

The above considerations also explain the much weaker 
gravitational wave signal of model AlBlGl in case of rel- 
ativistic gravity (see upper panel in Fig. H). In general 
relativity the central parts of the inner core collapse to 
almost 2.5pnuc, and rapidly settle down to an equilibrium 
state with constant density after a single bounce (middle 
panel). Since the matter is more concentrated towards the 
center, the density is relatively small in the outer parts of 
the core (lower panel). The fluid motion in the outer re- 
gions of the core near r = 13.0 km, where the shock forms, 
is not in phase with the fluid motion in the inner regions 
(see two lower panels in Fig. In Newtonian gravity 
the inner core is much less compact. It has a lower cen- 
tral density (middle panel), but a much higher density at 
r — 13.0 km (lower panel). In this case the entire inner 
core oscillates coherently with large amplitude motions, 
which give rise to a large time variation of the quadrupole 
moment and thus a large gravitational wave signal, the 
maximum amplitude being almost twice as large as the 
corresponding one in relativistic gravity (see Table |^). 

3.3. Change of the collapse dynamics 

In the comprehensive study of Newtonian rotational core 
collapse carried out by Zwerger & Miiller (1997), mod- 



els with multiple bounces were observed quite often. 
Relativistic gravity can have a qualitative impact on the 
dynamics of these models. If the density increase due to 
the deeper relativistic potential is sufflciently large, a col- 
lapse which is stopped by centrifugal forces at subnu- 
clear densities (and thus undergoes multiple bounces) in 
a Newtonian simulation, becomes a regular, single bounce 
collapse in relativistic gravity. 

Only two of the 26 models of our sample, A2B4G1 
and A3B4G2, show unambiguous multiple bounces in rel- 
ativistic gravity. These two models share some common 
features. In the relativistic simulations their peak cen- 
tral densities associated with the individual bounces are 
much larger (up to a factor 7) compared to those of the 
Newtonian runs (Table |l|). Consequently, the time elapsed 
between subsequent bounces decreases by a factor of up 
to 4. These type II models also show the effect of density 



crossing (see Section 3.2). As a result, despite the increase 
of the central density and the decrease of the oscillation 
timescale (giving rise to a larger average acceleration of 
the inner parts of the core), the maximum gravitational 
wave amplitude is lower for the relativistic models. 
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Fig. 9. Effects of a coherent motion of the inner core 
on the gravitational wave signal in the relativistic (solid 
lines) and Newtonian (dashed lines) simulation of model 
AlBlGl. In the upper panel the evolution of the gravi- 
tational wave amplitude A^q is plotted. The middle and 
lower panels show the evolution of the central density pc 
and of the equatorial density prc evaluated at a radius 
r = 13.0 km, respectively. Contrary to the relativistic sim- 
ulation, the inner core shows large amplitude oscillations 
in Newtonian gravity, which causes a much larger gravita- 
tional wave signal. The vertical dotted lines indicate the 
time of bounce and the horizontal dotted line in the 
middle panel marks nuclear matter density Pnuc- 



When the central densities become very large, rela- 
tivistic gravity can even change the collapse dynamics 
qualitatively, i.e. the collapse type is altered with respect 
to the Newtonian case (see Table |]). We have observed 
mostly transitions from a Newtonian type II (multiple 
bounce collapse) to a relativistic type I (regular collapse) 
case, and less frequently transitions from or to type I/II. 
A typical example is model A1B3G1 (Fig. |l^). The change 
of type is clearly visible both in the evolution of the cen- 
tral density (upper panel) and of the gravitational wave 
signal (lower panel). 

Centrifugal forces cause most of the type II multiple 



bounce models of Zwerger & Miiller (1997) to bounce 
at densities below pnuc- Only a few models collapse to 
nuclear matter density, the largest maximum density at 
bounce being Pmaxb — 3.09 x 10^^ g cm^'^ in the case of 
model A4B1G2. Accordingly, in the two relativistic type II 
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Fig. 10. Evolution of the central density pc (upper panel) 
and the gravitational wave signal amplitude A^q (lower 
panel) in the relativistic (solid lines) and Newtonian 
(dashed lines) simulation of model A1B3G1. The verti- 
cal dotted lines mark the times of peak central density for 
each bounce, and the horizontal dotted line in the density 
plot marks nuclear matter density Pnuc- 



models the maximum bounce density does not exceed this 
threshold. Whenever the deeper relativistic gravitational 
potential drives the maximum density in the core beyond 
the threshold of pmaxb ^ 3 x 10^'* g cm~'^, a type transi- 
tion occurs (see Table ||) . 

A collapse type transition does not only alter the sig- 
nal waveform qualitatively, but also has important conse- 
quences for the maximum gravitational wave signal ampli- 
tude, as e.g. in the case of model AlBlGl (Fig. ^). In this 
Newtonian type I/II model the concurrence of both high 
densities and a coherent motion of the inner core yields a 
large gravitational wave amplitude of ~ 1700 cm. As the 
stronger relativistic gravitational potential leads to the 
almost immediate formation of a new equilibrium state, 
efficiently suppressing coherent large amplitude motions 
of the inner core, the maximum signal amplitude is only 
943 cm, which is 44% less than in the Newtonian case. 

In Newtonian gravity and for an adiabatic index 71 — 
1.325 or 71 = 1.320 (i.e. close to the initial value 7ini — 
4/3) even models with rather small rotation rates /3rotini 
and moderately differential rotation (i.e. larger values of 
A) show multiple bounces. Our simulations demonstrate 
that relativistic gravity drastically narrows the region 
in the parameter space of /^rotini and A where multiple 
bounce models are possible. Fig. |l^ shows that only mod- 
els with high rotation rates and moderate differential ro- 
tation undergo an unambiguous type II multiple bounce 
collapse. 
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Fig. 11. Parameter space spanned by the rotation law 
parameter A and the rotation rate /Srot for multiple 
bounce collapse models in relativistic (normal font) and 
Newtonian (italic font) simulations. Filled circles mark 
models which are of multiple bounce type both in rela- 
tivistic and Newtonian gravity, and empty circles mark 
those models which show multiple bounces in Newtonian 
gravity only. 



3.4. Rapidly and highly differentially rotating core 
collapse 

When rapidly (large /Jrotini) and strongly differentially 
(small A) rotating initial models collapse on a short 
timescale due to a low adiabatic exponent 71 , angular mo- 
mentum conservation results in a considerable spin-up of 
the core. Consequently, the increasing centrifugal forces 
stop the collapse before or slightly above nuclear matter 
density, as in type II multiple bounce models. However, 
the maximum density is reached off-center, the topology of 
the density distribution in the core being torus- like. Some 
initial models (e.g. A4B4 and A4B5) already exhibit a 
toroidal structure. During collapse the toroidal character 
of such a model is significantly enhanced, and the max- 
imum density pmax can exceed the central density pc by 
up to two orders of magnitude. Other rapidly and highly 
differentially rotating initial models, like A3B5, A4B1 and 
A4B2, in general do not change their initial oblate density 
stratification during the evolution. However, for particular 
values of the parameter 71, they also develop a toroidal 
topology during the infall phase. 

Due to their deeper gravitational potential, relativistic 
configurations tend to be more compact compared to the 
corresponding Newtonian ones (see Section 3.2). This also 
holds for relativistic models with a toroidal density strat- 
ification, i.e. for models rotating rapidly and strongly dif- 
ferentially. These models have off-center density maxima 
which exceed those of their Newtonian counterparts (see 
Table 0). In some borderline models which rotate moder- 
ately fast and not too differentially, the core develops a 
toroidal density stratification only in Newtonian gravity 
and possesses an oblate spheroidal density stratification 
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Fig. 12. Equatorial density profiles of tlie moderately 
rapidly rotating model A3B3G5 in case of relativistic 
(solid line) and Newtonian (dashed line) gravity, depicted 
after ring-down. In relativistic gravity tfie core has an 
oblate density structure. The less deep Newtonian gravi- 
tational potential leads to the development of a toroidal 
density configuration. The horizontal dotted line marks 
nuclear matter density pnuc- 



in relativistic gravity. A prototype of such a case is model 
A3B3G5 (see Fig. 

In both relativistic and Newtonian simulations of 
rapidly and highly differentially rotating cores, the grav- 
itational wave amplitudes are generally large, because 
of the large quadrupole moments of toroidal configura- 
tions (see Table |l|). If the rotation rate is modest (models 
A4B1 and A4B2), the density in the torus at bounce can 
reach high supranuclear densities in relativistic gravity. 
The characteristic frequencies of the gravitational wave 
signal of these toroidal models become particularly high, 
as they change to regular collapse type models in relativis- 
tic gravity. The combination of a high density torus and 
large accelerations (reflected by the higher frequencies) 
yields very large wave amplitudes up to almost 3400 cm, 
which is significantly larger than in Newtonian gravity. 
However, in the case of the extremely rapidly rotating 
models A4B4 and A4B5, centrifugal forces surmount rel- 
ativistic gravity. The collapse type does not change, and 
the bounce is partially due to centrifugal forces (type I/II 
collapse). Consequently, the densities at bounce are less 
extreme and the average accelerations are smaller than in 
less rapidly rotating toroidal models. Therefore, the max- 
imum signal amplitudes jAfg jmax of models A4B4 and 
A4B5 are smaller than in the corresponding Newtonian 
models. 

Compared to models with a moderate amount of ro- 
tation, the propagation of the shock wave is drastically 
different in a rapidly and differentially rotating model, as 
e.g. in model A4B5G5 (Fig. p^ ). In this model the den- 
sity maximum is located off-center already during collapse 
(left panel). After core bounce the inner core settles down 
to an approximate equilibrium state with a distinctive 
toroidal density stratification surrounded by an equatorial 



disk from which matter is accreted onto the torus (middle 
panel). In this accretion disk the shock encounters rela- 
tively dense matter with large infall velocities between 0.1c 
and 0.15c, and thus it propagates outward comparatively 
slowly along this direction. On the other hand, the region 
close to the rotation axis has been emptied of matter dur- 
ing the collapse, the average density being orders of mag- 
nitude lower than close to the equatorial plane. Therefore, 
here the shock propagates out very rapidly (right panel). 
As a result, the flow velocity in the meridional plane ex- 
hibits a jet-like structure, with matter being ejected from 
the inner core predominantly along the rotation axis. 

3.5. Evolution of the rotation rate 

The influence of relativistic gravity on the evolution of 
the rotation rate during core collapse is demonstrated in 
Fig. |l^, which shows the radial profiles of = ^/vsv^ at 
the time of core bounce for several representative mod- 
els. Both the maximum and average rotation velocities 
are higher in the relativistic models, and the maximum of 
ti0 is shifted to smaller radii. Note that this maximum is 
located near the edge of the compact remnant. The dif- 
ferences are most pronounced for model A1B3G1 (lower 
panel), which changes from multiple bounce collap se t o 
regular collapse due to relativistic effects (see Section 3.3). 
Model A1B3G1 also shows the largest rotation velocities 
of all our models, approaching 0.2c at r « 15 km dur- 
ing core bounce. This corresponds to a rotation period of 
T ^ 1.6 ms at the edge of the inner core. When collaps- 
ing from the same initial model (A1B3), but imposing a 
larger initial reduction of the adiabatic index to 71 = 1.280 
(i.e. model A1B3G5; upper panel in Fig. |lj), the spin-up 
is less and < 0.08c. This result can be explained as 
follows: In model A1B3G1 the reduced adiabatic index 
71 — 1.325 is very close to 4/3, i.e. a large fraction of the 
core collapses nearly homologously ( poldreich fc Webei 



1980). However, this does not hold for model A1B3G5, 



where the outer parts of the core fall much slower than in 
model A1B3G1, and hence also fall less far in radius un- 
til the time of bounce. Consequently, the spin-up of mass 
zones located initially at the same radial position in both 
models is less in model A1B3G5 than in model A1B3G1 
except for the very center of the core. 

Assuming angular momentum conservation for the 
core, the rotation rate /3rot scales like the inverse of t he 
radius of the core (see, e.g., Shapiro fc Tcukolsky (1983 )), 



Prot — 



Ere 



J2 



1 



■pot I 



APR i?' 



(2) 



where E^ot, ^'poti M, R and J are the rotational energy, 
potential energy, mass, radius and angular momentum of 
the core, respectively. 

Knowing the behavior of the evolution of (3rot during 
core collapse is important for several reasons. Firstly, the 
rotation state of the inner core after bounce determines 
the rotation rate and profile of the compact remnant pro- 
vided no significant fall-back of matter occurs later in 
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Fig. 13. Formation of a torus and shock-propagation in the very rapidly and highly differentially rotating model 
A4B5G5. The three snapshots show color coded contour plots of the logarithm of the density, log p„uc, together with 
the meridional flow field during the infall phase a,t t — 25.0 ms (left panel), shortly after the centrifugal bounce at 
t = 30.0 ms (middle panel), and at t = 35.0 ms (right panel) when the torus is surrounded by large scale flow vortices. 
The length of the velocity vectors is scaled to the maximum value in the plotted region. Note that the size of the 
displayed region and the color coding (which is given above each panel) vary from plot to plot. 



the evolution. Pulsar observations show that neutron stars 
can spin with ve ry high rotatio n periods of the order of 
milliseconds (see Lorimer (2001 ) and references therein). 
However, such millisecond pulsars are believed to be born 
not as fast rotators, but are thought to be old objects 
which are spun up via mass transfer from a companion star 



much later in their evolution. Observations of pulsars in 
young supernova remnants (see, e.g. Kaspi (200C| )) provide 
only upper limits on the rotational state of newborn neu- 
tron stars, because significant angular momentum losses 
may have occurred since their formation. Presently, the 
fastest known young pulsar is the 16 ms X-ray pulsar in 
the Crab-like supernova remnant N157B ( Marshall et al 
1998|)~^econdly, at sufficiently high rotation rates the core 
will become unstable against triaxial perturbations due 
to various mechanisms (see Btergioulas (1998| ) and refer- 
ences therein). Such nonaxisymmetric rotating configura- 
tions are an additional source of gravitational radiation. 

For uniformly rotating, constant density MacLaurin 
spheroids in Newtonian gravity, the threshold values for 
the development of triaxial instabilities on secular and 
dynamic timescales are /3rots cc ^ 13.8% a nd /3rotdyn ~ 
27.4%, respectively (see, e.g., [Tassoul (19"78|)). 



tion rates, i.e. for values /3rotscc > 0.14, when the configu- 
rations become more compact. 

The dynamic instability, on the other hand, is driven 
by hydrodynamics and gravity. The onset of the dynamic 
bar-mode instability has been investigated in full gen- 
eral relativ ity only recently ( Shibata et al. 200C ; ^aijo 
et al. 2001), because it requires fully nonlinear multidi- 



The secular instability, which is driven for example by 
gravitational radiation or viscosity, has been analyzed in 
compressible stars, both uniformly and differentially ro- 
tating, within linear perturbation theory in full general 



relativity by 


Bonazzola et al. (1996 




199S), Stergioulas 


& Frie. 


Iman (1998 


), and 


Yoshida & Eriguchi (199£). The 



analysis shows that /3rotsoc depends on the compactness 
of the star, on the rotation law, and on the dissipa- 
tive mechanism. Whereas /3rotsoc < 0.14 for the gravita- 
tional radiation-driven instability (for extremely compact 
or strongly differentially rotating stars /3rotsoc < 0.1), sec- 
ular instabilities driven by viscosity occur for larger rota- 



mensional simulations (for Newtonian investigations see, 
e.g., New et al. (2000), and references therein). The simu- 
lations show that in rclativistic gravity the critical rotation 
rates for the onset of the dynamic bar-mode instability are 
only slightly smaller (/3rotdyn ^ 0.24 — 0.25) than in the 
Newtonian case, and depend only very weakly on the de- 
gree of differential rotation (within the moderate range 
surveyed by the investigations). 

What do these general considerations concerning tri- 
axial instabilities imply for our models? Since the rota- 
tion rates are larger in rclativistic gravity than in the 
Newtonian case for the same initial conditions, we ex- 
pect that the criteria for the development of dynamic and 
secular triaxial instabilities are fulfilled in more models 
and for longer time intervals in rclativistic gravity. To see 
whether this statement is correct we consider the evolution 
of the rotation rate /3rot for models A2B4G1, A4B4G4 and 
A4B5G4 (see Fig. |l^). In all these models Prot reaches its 
maximum value at the time of bounce. If the core experi- 
ences multiple bounces, this is reflected in multiple peaks 
of Prot (see upper panel). As expected, the rotation rate in 
the rclativistic models exceeds that of the corresponding 
Newtonian ones. This holds both for the maximum value 
and the average value of the rotation rate. 

Using the Newtonian thresholds /Jrotsoc and /3rotdyn 
(see discussion above) to judge whether a model is prone 
to developing triaxial instabilities, we see that rclativistic 
effects do not cause an important change in the case of 
model A2B4G1 (upper panel of Fig. |l^, because /3rot > 
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Fig. 14. Radial profiles of tlie rotation velocity in 
tfie equatorial plane plotted shortly after core bounce 
obtained in the relativistic (solid lines) and Newtonian 
(dashed lines) simulation of models A1B3G5 (upper 
panel), A2B4G1 (middle panel) and A1B3G1 (lower 
panel). The vertical dotted lines mark the radius of max- 
imum rotation velocity. 
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Fig. 15. Evolution of the rotation rate /3rot in the relativis- 
tic (solid lines) and Newtonian (dashed lines) simulations 
of models A2B4G1 (upper panel), A4B4G4 (middle panel) 
and A4B5G4 (lower panel). The horizontal lines mark the 
critical rotation rates above which MacLaurin spheroids 
are secularly {Prot sec) and dynamically P^ot dyn unstable 
against triaxial perturbations in Newtonian gravity. 



Aotsoc for less than 10 ms. On the other hand, the be- 
haviour of the rapidly rotating model A4B4G4 (middle 
panel) is quite different. Its rotation rate exceeds /3rotsoc 
near bounce time and remains above that value after core 
bounce in relativistic gravity (even exceeding /3rotdyn for 
a short time), while in Newtonian gravity it falls below 
Pi-ot sec again after ring-down. Since relativistic effects tend 
to lower the threshold value /3rot sec in case of the gravita- 
tional radiation-driven instability (see discussion above), 
this model demonstrates that Newtonian simulations can 
yield incorrect predictions about the development of sec- 
ular instabilities. On the other hand, in the extremely 
rapidly and differentially rotating model A4B5G4 (lower 
panel), secular instabilities are likely to develop in both 
Newtonian and relativistic gravity. However, instabilities 
on a dynamic timescale should only appear in the latter 
case, where /3rot > /3rotdyn for almost 4 ms. We note that 
for a neutron star with an average density of ~ pnuc and 



a radius of ^ 20 km, the dynamical timescale (= sound 
crossing timescale) is roughly 1 ms. 

Note that in this discussion /3rot is the ratio of rota- 
tional energy to potential energy of the entire core. Strictly 
speaking, the criteria for triaxial instabilities should only 
be evaluated for the rotation rate /3rot ic of the inner core 
(which eventually becomes the compact remnant), be- 
cause the outer core is a separate entity as far as the 
dynamics is concerned. However, due to the rapid drop 
of density and rotation velocity observed in all models at 
the outer boundary of the inner core, the outer core con- 
tributes only little to the total rotation rate. Thus, we can 
safely assume that /3iotic ~ Prot- This behavior has also 
be en observed in the Newtonian core collapse simulations 
of Zwerger (1995 ). 

Despite the knowledge of the evolution of /3rot during 
core collapse, predictions about the actual development 
of triaxial instabilities are difficult. The thresholds orig- 
inate from a linear stability analysis and involve restric- 
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tions (e.g. homogeneous density, uniform rotation, etc.). 
Unfortunately, we cannot simulate the actual nonlinear 
growth of triaxial instabilities, because our code is re- 
stricted to axisymmetric flows. Therefore, we have begun 
to extend our code to three spatial dimensions. The results 
of the 3D simulations planned to be performed with the 
extended code will be presented in a future publication. 
Still, at present our results may help to predict the be- 
havior and consequences of /?iot of fully three-dimensional 
models ( [Houser et al. 1994 ; Smith et al. 1996; Rampp et al 
1998()7] 

For large values of /3rot the neutron star is highly non- 
spherical. In this regime we expect deviations of the CF 
metric from the exact metric in the range of a few percent 
(see Section 5.5 in Paper I). Additionally, in the equation 
for /3rot (see Sections 1 and 5.4 in Paper I), the poten- 
tial energy Epot is determined by the (small) difference 
between the gravitational mass -Egiav and proper mass 
-^'proper- Evcn for strongly gravitating systems like neutron 
stars, this difference is only about 10% of the individual 
terms. Thus, if the values for -Egrav and improper are subject 
to an error of a few percent in the CFC approximation, 
we expect (3rot to deviate by several 10% relative to its 
value in an exact spacetime. In simulations of rapidly ro- 
tating neutron stars in equilibrium we could observe rel- 
ative deviations between /3rot exact and AotCFC of about 
15%. Owing to the uncertainties with the stability thresh- 
olds and the calculation of the rotation rate in the CFC 
approximation, the above predictions from the evolution 
of /3rot for the actual development of triaxial instabilities 
should be regarded as an estimate. 



4. Gravitational wave emission 

4.1. Gravitational wave energy spectrum and emission 

A detailed discussion of the spectra of gravitational waves 
emitted during Newtonian rotational core collapse, and 
their dependence on the collapse type and dynam ics can 
be found in the work of Zwerger fc Miiller (1997 ). In the 
quadrupole approximation and in axisymmetry, the grav- 
itational radiation field is solely determined by the am- 
plitude Afp (see Appendix The total energy radiated 
away by gravitational waves is then given by 



rad tot 



327r 



dt 



dt. 



(3) 



By replacing the gravitational wave amplitude in the time 
domain Afp (i) by its Fourier transform in the frequency 
domain Afp (i/), the total radiated energy i^radtot can be 
expressed as a frequency integral: 
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The spectral energy distribution is then 
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Fig. 16. Spectral energy distribution di?rad/di/ of the 
gravitational wave signal in the relativistic (solid lines) 
and Newtonian (dashed lines) simulations of models 
A1B3G3 (upper panel), A3B3G1 (middle panel) and 
A1B3G1 (lower panel). The energy spectrum of model 
A1B3G3, with a strong low frequency part and a weaker 
high frequency part, is typical for the signal of a regu- 
lar collapse. Model A3B3G1 is a multiple bounce collapse 
model in both Newtonian and relativistic gravity. In its 
energy spectrum a broad low frequency peak dominates. 
In model A1B3G1 relativistic effects change the collapse 
type from multiple bounce collapse to regular collapse. 
The energy spectrum reflects this change. The left (right) 
vertical dotted lines mark the maxima of the energy spec- 
trum in Newtonian (relativistic) gravity. Due to higher av- 
erage central densities in the relativistic simulations, the 
maximum of the energy spectrum is shifted to higher fre- 
quencies. 



As a characteristic spectral measure of the gravita- 
tional wave signal we define the frequency t'max at which 
the energy spectrum d£^r ad/di^ has its maximum (see also 
[Zwerger fc Miiller (199?!) ). 



As discussed in Section 3.2, the average and peak cen 



tral densities of our relativistic models are larger than 
those of the Newtonian ones. This leads to an increase 
of the oscillation frequencies of the inner core during the 
ring-down phase in a regular collapse situation, and to a 
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decrease of the time interval between consecutive bounces 
in multiple bounce models. Therefore, the gravitational 
wave spectrum is shifted towards higher frequencies in 
the relativistic models. This shift can be as large as a 
factor 5 in frequency (see Table |^ and lower panel of 
Fig. |6|). The models with the largest frequency shift are 
those where the deeper relativistic gravitational potential 
changes a multiple bounce collapse into a regular collapse 
(see Section 3.3). We find that in Newtonian gravity only 
three models have a higher characteristic frequency z^max 
than the corresponding relativistic ones. However, the en- 
ergy spectra of these models show many local maxima 
near z^max having almost the same spectral energy den- 
sity. Fitting the energy spectra of the three models by a 
smooth function yields again a higher value of Vniax in case 
of relativistic gravity. 

The energy spectra of three models are plotted in 
Fig. 16 . They are obtained by Fourier transforming the 



4.2. Gravitational wave amplitudes and frequencies 

The distribution of the dimensionless signal amplitude 
seen by an observer located in the equatorial plane of the 
source, i.e. the gravitational wave strain (see Eq. (A. 3) 
with 9 = 7r/2) 
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(6) 



is shown in Fig. for all 26 models, scaled to a distance 
of r = 10 kpc from the source. The corresponding wave 
amplitudes are listed in Table |l]. 



gr avita ltional signal. A Welch filter is applied in the Fourier 
transformation to reduce the noise caused by the finite 
length of the wave signal ( Press ct al. 1992 ). To compen- 
sate for the power loss due to filtering, the spectrum is 
rescaled by the requirement that the total energy in the 
time domain, Eq. (^ equals that of the frequency do- 
main, Eq. (||). The energy spectra are similar for both 
relativistic and Newtonian simulations in cases where the 
collapse type does not change (see upper and middle panel 
in Fig. |l^). On the other hand, when a Newtonian model 
exhibiting multiple bounces changes to a regular collapse 
model in relativistic gravity, like e.g. model A3B3G1, the 
shape of its energy spectrum is shifted significantly (see 
lower panel in Fig. |l6|). 

By integrating the spectral energy distribution, 
Eq. (|^), over all frequencies, one obtains the total energy 
radiated away by the system in gravitational waves, as 
given in Eq. (^). However, when analyzing numerical data 
this method is inappropriate, because high frequency noise 
contributes significantly to i?iadtot- For most models it is 
impossible to unambiguously define a cutoff frequency be- 
yond which there exists only noise, as -Erad does not drop 
sharply at some specific frequency. On the other hand, ap- 
plying a filter (like the Welch filter used in Fig. ^6|) causes 
a leakage of spectral energy density, and thus leads to a 
underestimate of -Eradtot- Therefore, we use Eq. (||) to cal- 
culate Ey-ad tot ■ 

We find that the total energy (measured in units of 
Mqc^) radiated in form of gravitational waves lies in the 
range 3.0x10"^° < -B,.adtot < 3.7 x 10""^ for the relativistic 
models, and in the range 3.8 x 10"^" < i?radtot < 1-5 x 
10~^ for the Newtonian ones. Averaged over all models, 
-Eradtot = 8.2 X 10^® in the relativistic case, and -Eradtot = 
3.6 X 10^^ in the Newtonian one, which corresponds to an 
increase of 128% due to relativistic effects. 



A comparison of Fig. |l^ with Fig. 8 in Zwerger & 



Mlillcr (1997) proves that our sample, encompassing 25 



of the 78 models considered in the parameter study of 
[Zwerger fc Miiller (1997 ) - we have simulated the addi- 
tional model A3B2G4soft - covers almost the same range 
of wave amphtudes (4 x lO^^i < /jTT < 3 ^ IQ-^O) and 
frequencies (60 Hz < < 1000 Hz). Thus, relativistic 
effects have no significant impact on the gross distribu- 
tion of expected signals in the amplitude-frequency dia- 
gram. Only in the high frequency, low amplitude region 
of the diagram, i.e. for models which collapse rapidly and 
have type I or HI waveforms, the weaker signal amplitudes 
cause a discernible change in the model distribution. The 
bulk of the models falls into the range of detectability 
of the LIGO I interferometer (Gustafson et al. 1999) pro- 
vided the source is located at a distance of less than 10 kpc 
(see Fig. |l^). 

For all models which are of the same type in both 
Newtonian and relativistic gravity, the gravitational wave 
signal has a lower amplitude. When selecting those 13 
models which are of type II in Newtonian gravity (see 
Section p. 1.2 ), the influence of relativistic effects on the 
signal amplitude is not as clearcut. The two relativis- 
tic models which remain of type II (see Table |l|) yield 
smaller maximum amplitudes I^^q jmax- If the collapse 
type changes due to the higher central densities reached 
in relativistic gravity, some models exhibit significantly 
higher signal amplitudes than in Newtonian gravity (e.g. 
-1-84% for model A1B3G1), whereas others possess smaller 
signal amplitudes (e.g. —22% for model A3B1G1). In mod- 
els which change their collapse type, the increase or de- 
crease of the signal strength, and thus the change of lo- 
cation in the amplitude-frequency diagram, depends very 
strongly on the rotational state of the initial model and 
on the value of 71. However, as in all other models, rel- 
ativistic gravity always causes a rise of the characteristic 
signal frequencies. 

The overall behaviour of the 13 models which are of 
type II in Newtonian gravity, can be summarized by con- 
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Fig. 17. Prospects of detection of the 
gravitational wave signal from axisym- 
metric rotational supernova core col- 
lapse in relativistic (black filled cir- 
cles) and Newtonian (red unfilled cir- 
cles) gravity. The figure gives the (di- 
mensionless) gravitational wave ampli- 
tude h^"^ and the frequency range for 
all 26 models. For a source at a dis- 
tance of 10 kpc the signals of all models 
are above the burst sensitivity of the 
LIGO I detector (except for some low 
amplitude, high frequency models), and 
well above that of the LIGO II interfer- 
ometer. The burst sensitivity giv es the 
r.m.s. noise amplitude /irms ~ \/ vS{v) 
over a bandwith of width ;/ at a fre- 
quency V for the instrument noise power 
spectral density Slv). The error bars 
mark the frequency range inside which 
the spectral energy density is within 
50% of its peak value. 
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Fig. 18. Impact of the frequency and 
wave amplitude shift on the sig- 
nal detectability for relativistic models 
which experience multiple bounces in 
Newtonian gravity. In the "best case", 
both amplitude and frequency increase 
in relativistic gravity (a). However, on 
average only the frequency increases 
(b). In neither case the shift is as steep 
as the high frequency sensitivity curve 
of an interferometer detector, and thus 
a relativistic model may fall out of the 
sensitivity window of the detector, if 
the Newtonian model is just within de- 
tection range. Note that the vertical 
scale is arbitrary. Symbols and error 
bars as in Fig. O. 



sidering the arithmetic mean of their peak amphtudes and 
peak frequencies: 
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Our simulations show that while the average gravitational 
wave signal amplitude K^"^ does not change significantly, 
the frequency more than doubles. 

This large frequency shift can have important con- 
sequences for the prospects of detection of gravitational 
waves. If the signal of a particular model is close to the 



detection threshold of a detector, as shown in Fig. [T^ 
for the 13 models which experience multiple bounces 
in Newtonian gravity, the influence of relativity can be 
twofold. In the "best case scenario" the signal is shifted to 
both higher frequencies and amplitudes (case a). However, 
more frequently the location of the model will only shift 
towards higher frequencies in the amplitude-frequency di- 
agram (case b - see Eq. (^), and therefore in a borderline 
case the signal may actually leave the sensitivity window 
of the detector. 



5. Summary and conclusions 

We have presented hydrodynamic simulations of relativis- 
tic rotational supernova core collapse in axisymmetry and 
have analyzed the gravitational radiation emitted by such 
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an event. We have simulated the evolution of 26 models 
in both Newtonian and relativistic gravity. Collapse of the 
initial configurations, which are differentially rotating rel- 
ativistic 4/3-polytropes in equilibrium, is induced by de- 
creasing the adiabatic index to some prescribed fixed value 
7i with 1.28 < 7i < 1.325. The stiffening of the equation 
of state at nuclear matter density and the thermal pressure 
in the matter heated by the prompt shock are simulated 
by means of a simplified equation of state consisting of a 
polytropic and a thermal part. Microphysics like electron 
captures and neutrino transport is neglected. 

Our simulations show that the three different types 
of rotational supernova core collapse and gravitational 
waveforms identified in previous Newtonian simulations 
(regular collapse, multiple bounce collapse, and rapid col- 
lapse) are also present in relativistic gravity. We find, how- 
ever, that rotational core collapse with multiple bounces, 
which occurs when the collapse is only or predominantly 
stopped by centrifugal forces, is possible only in a much 
narrower parameter range in relativistic simulations. Only 
two of the 26 models of our sample show unambiguous 
multiple bounces in relativistic gravity. The peak central 
densities associated with the individual bounces in such 
models are much larger (up to a factor of 8) in the rela- 
tivistic models compared to the corresponding Newtonian 
ones. Consequently, the time elapsed between subsequent 
bounces decreases by a factor of up to 4. 

We further find that relativistic gravity can have a 
qualitative impact on the dynamics. If the density increase 
due to the deeper relativistic potential is sufficiently large, 
a collapse which is stopped by centrifugal forces at subnu- 
clear densities (and thus undergoes multiple bounces) in 
a Newtonian simulation, becomes a regular, single bounce 
collapse in relativistic gravity. A collapse type transition 
also has important consequences for the maximum gravi- 
tational wave signal amplitude. 

For all relativistic models the peak density at bounce 
is larger than in the corresponding Newtonian models, the 
relative increase reaching ~ 700% in some cases. This also 
holds for the maximum density of the compact remnant af- 
ter ring-down in case of instant formation of a stable equi- 
librium state (types I and III). Nevertheless, we find that 
only six models have also larger maximum signal ampli- 
tudes than their Newtonian counterparts. The maximum 
signal amplitudes of all other relativistic models arc (up to 
57%) smaller than those of the corresponding Newtonian 
ones. The reduced maximum signal strength can be ex- 
plained by the fact that the quadrupole amplitude is de- 
termined by the bulk motion of the core rather than just 
by the motion of the densest mass shells. Therefore, a core 
which is more condensed in the center can give rise to a 
smaller gravitational wave signal than a core which is less 
centrally condensed, but which is denser and moves faster 
in its outer regions. However, the fact that all relativistic 
models show an increased central and reduced peripheral 
density of the collapsed core relative to their Newtonian 
counterparts does not necessarily imply a weaker gravi- 
tational wave signal, as in 6 of 26 relativistic models the 



maximum gravitational wave amplitudes are larger than 
in the corresponding Newtonian models. In all these six 
cases the collapse changes from type II (or type I/II) to 
type I when changing from Newtonian to relativistic grav- 
ity. 

Overall, the relativistic models cover almost the same 
range of gravitational wave amplitudes (4 x 10^^^ < 
/I'^T < 3 X 10-20) and frequencies (60 Rz < v < 1000 Hz) 
as the corresponding Newtonian ones. For all models 
which are of the same collapse type in both Newtonian and 
relativistic gravity, the gravitational wave signal is of lower 
amplitude. If the collapse type changes, either weaker or 
stronger signals are found in the relativistic case. Averaged 
over all models, the total energy radiated in form of grav- 
itational waves is 8.2 x 10^^ Mqc^ in the relativistic case, 
and 3.6 x 10~^ Mqc^ in the Newtonian one. 

For a given model, relativistic gravity can cause a large 
increase of the characteristic signal frequency of up to a 
factor of five. This frequency shift is an effect common to 
all models, and can have important consequences for the 
prospects of detection of gravitational waves. If the signal 
of a particular model is close to the high-frequency detec- 
tion threshold of a detector, the influence of relativity can 
be twofold. In the "best case scenario" both the signal fre- 
quency and amplitude increase; however, more often only 
the signal frequency increases. In either case, in a border- 
line situation the signal may actually leave the sensitivity 
window of the detector, as the high-frequency detection 
threshold of a laser interferometer rises steeply with in- 
creasing frequency. Hence, within the limits of our phys- 
ical model of rotational core collapse, relativistic effects 
on average tend to decrease the prospects of detectabil- 
ity of gravitational radiation by interferometric detectors. 
Nevertheless, the gravitational wave signals obtained in 
our study are within the sensitivity range of the first gen- 
eration laser interferometer detectors if the source is lo- 
cated within the Local Group. 

In our study we have restricted ourselves to simpli- 
fied supernova core models. The complicated microphysics 
(electron captures and other weak interactions, equation 
of state) and neutrino transport processes have all been 
absorbed in one model parameter (71). Although this is an 
extreme simplification, we nevertheless think that it does 
not change the qualitative nature of our results. In partic- 
ular, we do not expect that more realistic (axisymmetric) 
core collapse models will produce significantly stronger 
gravitational wave signals. Larger initial rotation rates 
and lower values of the (effective) adiabatic index 71, 
which could lead to stronger gravitational wave signals, 
are most likely excluded. This is also supported by the fact 
that the signal strengths of the microscopically more real- 
istic, but Newtonian models of Monchmeyer et al. (1991) 
agree well with those of our simplified models. To check 
whether these expectations are indeed correct we plan to 
incorporate a realistic equation of state and a simplified 
treatment of neutrino transport into our relativistic code 
in the near future. 
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In several models the rotation rate of the compact rem- 
nant exceeds the critical value where MacLaurin spheroids 
become secularly, and in some cases even dynamically 
unstable against triaxial perturbations (/3rotscc ~ 13.8% 
and /?rotdyn ~ 27.4%). It has been have pointed out that 
if such instabilities indeed occur, a much stronger grav- 
itational wave signal might be produced than by cores 
which remain axisymmetric (for a review, see e.g. Thorne 
(1995|)]. According to this idea, the core will be trans- 
formed into a bar-like configuration that spins end-over- 
end like an American football. One has further speculated, 
whether the core might even break up into two or more 
massive pieces, if /3rot > /Srotdyn- In that case the result- 
ing gravitational radiation could be almos t as strong a s 
that from coalescing neutron star binaries ( Thorne 1995| ). 
The strength of the gravitational signal sensitively de- 
pends on what fraction of the angular momentum of the 
non-axisymmetric core goes into gravitational waves, and 
what fraction into hydrodynamic waves. These sound and 
shock waves are produced as the bar or lumps, acting like a 
twirling-stick, plow through the surrounding matter. In or- 
der to investigate these effects, three-dimensional relativis- 
tic simulations of rotational core collapse have to be per- 
formed, which can follow the growth of non-axisymmetric 
instabilities (bar modes and gravitational wave driven un- 
stable inertial modes). To this end we are currently ex- 
tending our code to three spatial dimensions. 
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Appendix A: Gravitational wave extraction 

In the appendix we summarize the technical aspects re- 
lated to the gravitational wave extraction and the differ- 
ent methods we have used. As stated in Paper I, the use 
of a CF metric implies the removal of the gravitational 



radiation degrees of freedom. Therefore, in order to com- 
pute the waveforms we apply the Newtonian quadrupole 
formula. 

In that approximation and for an axisymmetric source 
in spherical coordinates, it can be shown (Thorne 1980) 
that the transverse traceless gravitational field has one 
independent component hj^ , which depends only on the 
quadrupole signal amplitude Afg , 
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spherical tensor harmonic T, 
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defined according to 
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the quadrupole radiation field hgg is given by 
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A.l. Different formulations of the quadrupole formula 



The amplitude Afg in Eq. ( A.3 ) is the second time deriva- 
tive of the mass quadrupole moment of the source. 
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where z = cos9, and k ~ IQtt^/-^ /^/T5. This formulation 
of the radiation field is known as the standard quadrupole 
formula. 

A direct numerical implementation of this formula 
is problematic, as the discretization of the second time 
derivative causes high frequency noise in the gravitational 
wave signal. The amplitude of this noise can be larger than 
the total amplitude of the signal. The problem is worsened 
by the fact that the density distribution in the integrand of 
Eq. ( |A.4| ) is weighted by r^, i.e. mass elements in the outer 
parts of the star - where the grid resolution is coarser - 
contribute strongly to the signal. 

To overcome this problem a standard practice is to 
replace the time derivatives by spatial derivatives (Finn 



1989; Blanchct ct al. 1990). This procedure leads to the 



stress formula for the gravitational wave amplitude: 

= k Jp(v^{3z^-l) + v^{2-3z^)-vl-6rVrVsVl^' 



-r^{3z^-l) + 3^Vl-z^]r^drdz, (A.5) 



dr 



89 



where $ is the Newtonian gravitational potential. This 
modified formulation of the quadrupole formula was used 
in our analysis. 

A.2. Ambiguity of the stress formula in the CFC 
approximation 

In a mathematically strict sense the quadrupole approxi- 



mation of the radiation field, Eq. (A.l), is only defined for 
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in Afg for rapidly rotating neutron stars in equilibrium, 
where the gravitational radiation should vanish. 

The Newtonian potential can also be computed ac- 
curately up to the first post-Newtonian order by setting 
gii ~ (1 — 2<i>) equal to the conformal factor This leads 
to 



(A.7) 
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The waveforms obtained from this expression also show 
an offset, but this time a negative one (see the dashed line 
in Fig. A.l). However, combining both expressions for <i>, 
Eq. ([A.6|) and Eq. (O), into 



Fig. A.l. Reduction of the offset in the gravitational wave 'i' weighted = 
amplitude A^q for model A3B2G4 with the weighting 
method applied to the stress formula. If the Newtonian 
potential $1 is used, the wave signal shows a large posi- 
tive offset (dashed-dotted line) . If the combination of met- 
ric components $2 is used, a large negative offset occurs 
(dashed line) . The weighted sum weighted of these two sig- 
nals yields a waveform with practically zero offset (solid 
line). This combined signal exhibits no high frequency 
noise, as its constituents have been obtained using the 
stress formula. 
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where a = const, is a parameter with an appropriate value, 
one is able to reduce the offset almost to zero in all core 



collapse models (see the solid line in Fig. A.l) 



a flat Minkowski spacetime where the gravitational waves 
propagate as a linear perturbation. However, the multipole 
amplitudes and the tensor harmonics, although defined in 
flat spacetime, contain information about the wave source, 
which is situated in a strong field and high velocity re- 
gion. In our simulations the quadrupole moment Afp is 
computed as a spatial integral over the relativistic matter 
configuration in a curved spacetime, as the metric f/^,y is 
not Minkowskian. Therefore, our methods to compute the 
radiation field are hampered by several ambiguities in the 
quadrupole formula. 

Whereas in the Newtonian formulation, which is based 
on Euclidean geometry, the definition of the radius coor- 
dinate r is no source of ambiguity, in the Riemannian ge- 
ometry of a curved spacetime the circumferential radius 
and the coordinate radius do not necessarily coincide. In 
our simulations we use the isotropic radial coordinate r in 
the computation of the quadrupole moment. Tests where 
we replaced r by the circumferential radius in the wave 
extraction produced almost identical waveforms. 

The spatial derivatives of the Newtonian potential 



present in Eq. ( A. 5 ) need particular attention. When sim- 



ply using Poisson's equation, 

A<I>i = 47rp, (A.6) 
which we solve by an expansion into Legendre polynomials 



(Zwerger 1995), the waveform shows an offset, particularly 



after core bounce (see the dashed-dotted line in Fig. A.l). 
For strong gravity, the Newtonian potential <I>i is therefore 
inappropriate for describing the radiation field. Another 
manifestation of this problem is a large constant offset 



